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Abstract 

Experimental  data  show  that  heterogeneous  nanostructure  of  solid  oxide  and  polymer  electrolyte  fuel  cells  could  be  approximated  as  an 
infinite  set  of  fiber-like  or  penny-shaped  inclusions  in  a  continuous  medium.  Inclusions  can  be  arranged  in  a  cluster  mode  and  regular  or 
random  order.  In  the  newly  proposed  theoretical  model  of  nanostructured  material,  the  most  attention  is  paid  to  the  small  aspect  ratio  of 
structural  elements  as  well  as  to  some  model  problems  of  electrostatics.  The  proposed  integral  equation  for  electric  potential  caused  by  the 
charge  distributed  over  the  single  circular  or  elliptic  cylindrical  conductor  of  finite  length,  as  a  single  unit  of  a  nanostructured  material,  has 
been  asymptotically  simplified  for  the  small  aspect  ratio  and  solved  numerically.  The  result  demonstrates  that  surface  density  changes  slightly 
in  the  middle  part  of  the  thin  domain  and  has  boundary  layers  localized  near  the  edges.  It  is  anticipated,  that  contribution  of  boundary  layer 
solution  to  the  surface  density  is  significant  and  cannot  be  governed  by  classic  equation  for  smooth  linear  charge.  The  role  of  the  cross-section 
shape  is  also  investigated.  Proposed  approach  is  sufficiently  simple,  robust  and  allows  extension  to  either  regular  or  irregular  system  of 
various  inclusions.  This  approach  can  be  used  for  the  development  of  the  system  of  conducting  inclusions,  which  are  commonly  present  in 
nanostructured  materials  used  for  solid  oxide  and  polymer  electrolyte  fuel  cell  (PEMFC)  materials. 
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1.  Introduction 

The  development  of  nanostructured  and  self-assembling 
materials  attracted  significant  attention  in  various  applica¬ 
tions,  such  as  energy  storage  and  energy  conversion  devices. 
However,  the  main  emphasis  of  nanoscale  studies  of  nano¬ 
materials  has  been  made  in  the  synthesis  of  these  materials, 
rather  than  in  the  area  of  a  rigorous  theoretical  approach. 
Thus,  some  of  the  existing  experimental  results,  such  as 
those  related  for  example  to  exceptionally  high  ionic  conduc¬ 
tivity  of  yttri a- stabilized  zirconia  ultrathin  films  [1]  cannot 
be  explained  on  the  basis  of  existing  mathematical  mod¬ 
els.  Mathematical  description  in  correlation  with  detailed 
structural  characterization  and  information  about  transport 
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properties  of  these  materials  could  open  new  opportunities 
and  significantly  extend  the  area  of  their  application. 

A  wide  variety  of  nanostructured  materials  are  commonly 
used  for  porous  catalyst  layers  of  solid  oxide  and  polymer 
electrolyte  fuel  cells  [2-6],  self-assembling  structures  such 
as  carbon  aerogels  [7],  carbon  nanotubes  [8],  gas  diffusion 
layers  [9,10],  and  gas  separation  membranes  [11,12].  These 
materials  can  be  assumed  to  be  a  combination  of  a  series 
of  conductive  inclusions  separated  by  a  non-conductive  gas 
medium. 

Various  mathematical  approaches  to  heterogeneous  struc¬ 
tures  have  been  suggested  [13-15].  One  of  these  applies  a 
straightforward  numerical  solution  of  a  3D  problem  by  means 
of  finite  element  or  boundary  element  (FE/BE)  code.  How¬ 
ever,  this  approach  is  feasible  only  in  the  case  of  a  single 
inclusion  or  a  finite  number  of  inclusions  and  is  not  robust 
when  for  a  structural  system  of  an  unlimited  number  of  inclu- 
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sions.  In  the  theory  of  elasticity  it  has  been  stated,  that  singular 
perturbation  of  a  domain  is  one  of  the  most  critical  issues  for 
solving  this  problem  numerically.  Thus,  a  direct  numerical 
technique  is  not  robust  when  we  consider  a  dynamic  process 
in  the  medium  with  a  singular  perturbation. 

Another  widely  used  approach  is  based  on  the  repre¬ 
sentative  volume  concept  when  “effective”  homogeneous 
parameters  of  an  heterogeneous  material  can  be  derived  from 
“representative”  volume,  i.e.  cubical  domains  containing  a 
finite  set  of  inclusions.  This  approach  to  a  heterogeneous 
material  allows  one  to  derive  constitutive  equations  of 
heterogeneous  medium  from  obtained  properties  of  the 
representative  volumes. 

The  goal  of  this  work  was  to  develop  a  simple  model  of  a 
thin  conducting  inclusion  in  a  dielectric  medium.  In  future, 
a  complete  theoretical  model  allowing  rigorous  predictions 
of  the  properties  of  nanostructured  materials  will  be  created 
on  the  basis  of  an  asymptotic  approach  to  a  single  conduc¬ 
tive  inclusion,  proposed  in  our  work.  In  addition,  it  should 
be  emphasized,  that  the  proposed  asymptotic  approach  can 
be  applied  to  a  dielectric  media  containing  systems  of  thin 
conducting  inclusions  with  various  positions,  orientations, 
dimensions,  configurations,  curvilinear  contours,  and  prop¬ 
erties. 


z 


Fig.  1.  Configuration  of  the  conductor. 


2.  Model  description 

In  order  to  solve  the  problem  of  a  single  cylindrical  inclu¬ 
sion,  an  asymptotic  simplification  is  used  with  the  assumption 
that  the  aspect  ratio  (5)  of  the  inclusion  is  small,  i.e.  two 
transversal  dimensions  are  much  smaller  than  the  length. 
The  cross-section  may  differ  from  the  circle;  the  electrostatic 
charge  is  distributed  over  the  surface  of  the  conductor.  To 
derive  an  equation  for  the  unknown  density,  the  original  inte¬ 
gral  equation  for  the  potential  is  asymptotically  simplified. 
The  original  3D  problem  is  decomposed  into  a  combination 
of  two  problems  of  reduced  dimension.  The  first  one  is  gov¬ 
erned  by  the  integral  equation  for  the  unknown  linear  charge 
(net  charge  in  the  cross-section)  over  the  axis  of  the  conduc¬ 
tor.  The  second  one  is  a  2D  electrostatic  problem  governed  by 
the  integral  equation  over  the  boundary  contour  of  the  cross- 
section.  The  derived  asymptotic  equations  for  the  unknown 
density  are  valid  at  the  points  located  far  away  from  the  edges. 

In  order  to  determine  the  charge  density  near  the  edge, 
the  boundary  layer  equation  is  derived.  Uniform  circular  and 
elliptic  cylinders  are  used  as  examples.  Results  for  a  circular 
cylinder  correlate  well  with  classic  results  in  the  middle  part 
of  the  domain;  contribution  of  the  boundary  layer  solution  is 
significant  near  the  ends;  influence  of  cross-section  shape  is 
investigated. 

2.7.  Formulation  of  the  electrostatic  problem 

The  perfect  conductor  with  a  total  electrostatic  charge  ( Q ) 
is  placed  into  an  infinite  dielectric  medium  with  dielectric 


constant  s.  Charge  is  assumed  to  be  distributed  over  the  sur¬ 
face  (S)  of  the  conductor.  Configuration  of  the  conductor  is 
shown  in  Fig.  1,  where  r  is  the  mid-line  contour  (closed 
or  unclosed),  y  the  boundary  contour  of  the  cross-section 
(contour  y  can  differ  from  the  circle).  The  unknown  surface 
density,  p,  of  the  charge  is  governed  by  the  classic  equation 
of  electrostatics: 


1  f  f  p(U ,  / ,  z')dS(x' ,  / ,  zQ 

4  ns  s  fix-  vf  +  (y-  y')2  +  (z-  z')2 


(2.1.1) 


Q  =  JJ  p(x',y',z)dS(x',y',z).  (2.1.2) 

Here  V  stands  for  the  electric  potential  of  the  conduc¬ 
tor,  v,  y,  z  the  Cartesian  coordinates,  and  Q  is  the  total 
charge. 

Asymptotic  analysis  of  the  integral  in  the  right-hand 
side  of  Eq.  (2.2.1)  has  been  conducted  in  potential  the¬ 
ory  (electrostatics/electrodynamics,  acoustics,  hydrodynam¬ 
ics,  theory  of  thin  aerofoil,  and  material  science/composites) 
[3].  In  this  paper,  we  apply  results  obtained  in  asymp¬ 
totic  analysis  of  singular  integrals  to  specific  problem  of 
electrostatics. 


2.2.  Axially  symmetric  problem 

We  consider  first  the  electrostatic  problem  for  a  circular 
cylinder  in  the  axially  symmetrical  electrostatic  field.  We  use 
an  integral  equation  (2.1.1)  and  assume  that  the  edge  area  is 
negligible  and  the  surface  S  means  the  lateral  surface  of  the 
conductor.  The  density  depends  upon  the  axial  coordinate  z, 
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p  =  p(z).  Eq.  (2.1.1)  can  be  rewritten  in  the  form: 


p(z')dz'dip 

\J(z  —  z')2  +  4 a2  sin2(<^/2) 


(2.2.1) 


ip  stands  for  polar  angle. 

The  unknown  surface  density  is  represented  as  a  sum  of 
the  “smooth”  part  describing  the  density  in  the  points  far 
away  from  the  ends,  and  two  boundary  layers  localized  near 
the  ends  z  =  L  and  z  =  —  L: 


p(z)  =  q  (f)  +  0 


L  —  z 


+  0 


L  +  z 


a 


a 


(2.2.2) 


2.2.7.  Equation  of  boundary  layer  near  the  edge 

After  solving  Eq.  (2.2.3)  for  the  smooth  part  q(s)  of  the 
density  we  can  calculate  density  corrected  for  the  bound¬ 
ary  layer  near  the  edge  z  =  L.  For  the  reason  of  sym¬ 
metry,  the  boundary  layer  near  the  opposite  edge  z  =  —L 
is  assumed  to  be  identical.  We  substitute  the  representa¬ 
tion  p(z)  =  q(z/L)  +  0((L  —  z)/a)  into  Eq.  (2.2.1)  under  the 
assumption  that  a  contribution  of  the  boundary  layer  near 
the  opposite  end  (z  =  —L)  is  negligible  near  the  end  z  =  L 
under  study.  Thus,  introducing  a  new  “short”  variable, 
r  =  (L  —  z)/a,  we  subtract  Eq.  (2.2.3)  for  a  smooth  part.  After 
some  asymptotic  replacements  and  limit  transition  5—^0, 
1/5  — >  oo  we  derive  the  boundary  layer  equation  in  the 
form: 


0(E) 


'2jT 


(t  —  r')2  +  4  sin2(<^/2) 


dr'  dip  =  —q(l) 


ip  sin  ip  dip 


0  \J x1  +  4  sin2 (ip/ 2)  x  +  \J  x1  +  4  sin2(<p/2) 


(2.2.5) 


The  smooth  part  is  associated  with  the  large  scale  L  of  a 
change  in  the  variable  z;  the  boundary  layers  correspond  with 
the  small  scale  a. 

Equation  for  the  smooth  part  q(s)  of  the  density  is 
derived  using  the  asymptotic  integral  equation  assuming  that 
a  <^L  —  |z|,  a  <^L,  8  =  a/L<Jf  1: 


2s y  _  il  qW)_  -  g(s)dj, 
a  J\  \s  —  s'\ 


+  q(s) 


In 


1  -  v 


47r/2 


X  JO 


ln(sin  ip)dip 


(2.2.3) 


2.3.  Elliptic  cylinder 

Consider  a  cylindrical  conductor  with  elliptic  cross- 
section  (ellipse  y).  Let  a  and  b  ( b<a )  stand  for  semi-axes 
of  the  ellipse.  Parameters  a  and  b  are  of  the  same  order  (the 
aspect  ratio  K  =  b/a  =  0(1)  and  both  are  much  smaller  than 
the  length  2L(a/L  =  o(l)).  Surface  density  p  depends  on  the 
longitudinal  coordinate  z  (associated  with  “long”  scale  L) 
and  polar  angle  ip  within  the  cross-section.  Contour  y,  polar 
angle  ip  and  Cartesian  coordinates  x,  y,  and  z  are  shown  in 
Fig.  2. 

Potential  equation  (2.1.1)  can  be  written  in  the  form: 


P(<P',  zf) 

\J(z-  z')2  +  r2(ip,  ip') 


7  /  /  f  A  7 

=  r  (ip,  ip)  =  a  (cos  ip  —  cos  (p  )  +  b Esin  ip  — 


’  /  v 

sin  ip  ) 


=  4  sin 


2  ip-ip 


a 


2  ■  2  (<P  +  <P'\  .  ,2  2  (V  +  V' 

sin  -  +  b  cos  - 

2  V  2 


(2.3.1) 


where  s-z/E  a  new  dimensionlesss  variable,  small  terms 
0(8  In  5)  are  neglected.  This  equation  is  to  be  solved  together 
with  equality  (2.1.2)  at  a  given  total  charge  Q.  The  derived 
asymptotic  equation  correlates  very  well  with  the  classic 
equation  of  Pocklington  [4]  for  a  thin  cylindrical  conductor: 

a  fL  p(zf)dzf 

V  —  —  /  ,  =.  (2.2.4) 

J-L  \/(z  -  z')2  +  a2 

These  equations  are  asymptotically  equivalent,  and  differ  in 
the  negligible  term  0(5  In  5).  Eqs.  (2.2.3)  and  (2.2.4)  allow 
one  to  calculate  “smooth”  density  at  the  points  far-away  from 
the  ends.  However,  smooth  density  may  contain  some  error 
near  the  ends  of  the  conductor. 


First,  we  look  for  smooth  surface  density  p  and  linear  charge 
P  (net  charge  in  the  cross-section)  as  functions  of  Cartesian 


Contour  y 


Fig.  2.  Contour  y,  polar  angle  <p,  and  Cartesian  coordinates  x,  y,  and  z- 
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coordinates  x,  y,  and  z: 


f  Z  IV 

a  a  LJ  ’ 


7i  y 


z 


a 


9  ’  7 

a  L 


/  b  .  , 

cos  cp  ,  -  sin  cp  , 
a 


dy(x',  y) 


sin  <p'+b2  cos  cpr  dcp' . 


(2.3.2) 


Further,  we  develop  combined  linear  charge  (7)  including 
boundary  layers  localized  near  the  ends  in  the  form: 


m  =  FlL)+*(Czl 


+  <z> 


L  +  z 


a 


(2.3.3) 


We  are  seeking  the  unknown  linear  charge  P(s )  distributed 
along  the  z-axis  (“smooth”  part  of  solution),  and  charge 
distribution  q  =  q(x/a ,  y/a ,  s)  over  the  contour  y  (of  the  cross- 
section  z  =  const).  Then  we  analyze  linear  charge  T(z)  as  a 
function  of  “short  scale”  coordinate. 


2.3.1.  Asymptotic  equation  for  the  smooth  density 

First,  we  substitute  this  representation  (2.3.2)  for  the  den¬ 
sity  to  Eq.  (2.3.1)  and  derive  the  asymptotic  of  potential 
V: 


AtXSV  = 


‘L  P(z'/L)  -  P(z/L )  ,  ,  .  „  fz  \  ,  4 (L2-z2) 


dz'  +  P  -  In 


L 


Z  Z 


L 


a " 


b 


q  cos  (p  ,  —  sin  (p 
a 


/x  2 


y  \  T  / 

2  /  2  / 

+  /c  (sin  <p  —  sin  <p  )  ]dy(yp  ). 


ln[(cos  <p  —  cos  cp  ) 


(2.3.4) 


The  structure  of  Eq.  (2.3.4)  gives  an  indication  that  we  have 
to  seek  the  density  q  in  the  degenerated  form: 


q  (cos  (p,  ^  sin  cp,  fj  =  4>((p)P  (|-)  ; 


>P(<p')dy(<p')  =  1. 


(2.3.5) 


Y 


Then  we  substitute  (2.3.5)  into  (2.3.4)  and  achieve  a  system 
of  four  equations  for  four  unknowns:  P(zJL ),  7/r( <p ),  V,  and  c. 


47r£F  = 


'L  P(z7L)  -  P(z/L)  ,  ,  .  „  /  Z  \  ,  4(L2— z2) 


L 


Z  Z 


dz'+P  (  ^  )  In 


/x2 


2  In  a-\-  /  ^Piyp  )  ln[(cos  (p  —  cos  (p  ) 


y 


?  2 
T /c  (sin (p  —  sin <p  )  ]d y(cp')  =  2 Inc; 


^(cpf)dy((p')  =  1; 


y 


: '  a 


dz'  =  Q 


(2.3.6) 


Thus,  separation  of  variables  is  completed.  The  first  equation 
of  the  derived  system  contains  information  on  3D  effects  of 


the  phenomenon  under  study.  The  second  equation  corre¬ 
sponds  with  the  2D  electrostatic  problem  of  charge  distri¬ 
bution  in  a  cylindrical  conductor  of  infinite  extent.  It  has  a 
shortcoming  of  being  a  2D  formulation;  2D  does  not  allow 
the  evaluation  of  3D  effects,  and  of  the  determination  of  the 
net  charge  P  in  the  cross-section.  Unknown  parameter  c  can 
be  considered  as  some  “effective”  radius  of  circular  cylinder 
with  the  same  linear  charge  P(z!L). 


2.3.2.  Boundary  layer  near  the  edge 

After  solving  Eq.  (2.3.6)  for  the  linear  charge  P(z/L) 
and  angular  distribution  ^( cp )  we  can  calculate  the  correc¬ 
tion  to  the  linear  charge  P  caused  by  the  boundary  layer 
near  the  edge  z  =  L.  This  boundary  layer  is  constructed 
under  the  assumption  of  validity  of  angular  distribution 
of  charge  (within  the  cross-section)  everywhere,  including 
at  the  vicinity  of  this  edge.  We  substitute  representation 
p((p,z)  =  [P(z/L)  +  <P((L  —  z)/a)]^P((p)  into  Eq.  (2.3.1).  (Con¬ 
tribution  of  the  boundary  layer  localized  near  the  opposite  end 
z  =  —  L  to  the  density  near  the  end  z  =  L  is  negligible.)  Then  we 
incorporate  the  new  “short”  variable  r  =  (L  —  z)/a  and  subtract 
the  first  equation  (2.3.7)  for  a  smooth  part  P(zJL )  from  Eq. 
(2.3.1).  After  some  asymptotic  replacements  and  limit  transi¬ 
tion  5^0;  1/5  — >  oo  we  obtain  the  boundary  layer  equation 
in  the  form: 


In  poo 

n<p')d y(<p')  / 
ro  Jo 


X 


0(Z')  ,  [2lZ 

!  ^=dr/=E(l)  / 

y  (r  —  x')2  +  (1  —  cos  cp')2  +  k2{ sin  cp)2 


x  ^((^ln 


2t 


r  +  \jx2  +  (1  —  cos  (p')2  +  k2{  sin  cp'f 


dy((p'). 


(2.3.7) 


3.  Results  and  discussion 

High-resolution  images  of  conductive  heterogeneous 
structures  are  presented  in  Figs.  3  and  4.  The  micrographs 
indicate  that  cylindrical  or  almost  cylindrical  fibers  are  ele¬ 
ments  of  the  structure;  the  fibers  within  the  structure  can  be 
assembled  in  three  typical  orders:  tree-like  order  (cluster), 
chaotic  order,  and  regular  order.  All  these  orders  can  be  gov¬ 
erned  by  asymptotic  equations  over  the  graph.  The  graph  can 
be  connected  or  contain  fragments  without  connection  with 
the  main  tree. 


3.1.  Circular  cylinder 


Equations  for  the  smooth  part  and  the  boundary  layer 
were  solved  numerically.  Calculated  dimensionless  density 
and  linear  charge: 


2j raLp(z) 

Q 


2jtaLp(z) 

ffsrC')  ds(x',  y,  z') 


(3.1.1) 
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Fig.  3.  Electrically  conducting  carbon  paper  used  as  a  gas  diffusion  layer  in 
PEMFCs;  inset  shows  the  same  material  at  higher  magnification. 


Fig.  4.  Electrically  conducting  carbon  cloth  used  as  a  gas  diffusion  layer  in 
PEMFCs;  inset  shows  the  same  material  at  higher  magnification. 


Fig.  6.  Dimensionless  density  calculated  from  asymptotic  equation  at  vari¬ 
ous  aspect  ratio  (a/L  =  0.05;  0.01;  0.001). 


is  presented  in  Figs.  5-7.  Fig.  5  presents  the  density  calcu¬ 
lated  from  the  equation  of  Pocklington  at  various  aspect  ratios 
(a/L  =  0.005;  0.001;  0.0001).  Fig.  6  presents  the  dimension¬ 
less  density  calculated  from  the  asymptotic  equation  (2.2.3) 
(without  boundary  layer)  at  various  aspect  ratios  (a/L  =  0.005 ; 
0.001;  0.0001).  Comparison  between  the  linear  charge  cal¬ 
culated  from  the  equation  of  Pocklington  (dotted  line),  Eq. 
(2.2.3)  (solid  line),  and  Eqs.  (2.2.3)  and  (2.2.5)  (dashed  line) 
at  a/L  =  0.05  is  presented  in  Fig.  7.  Calculations  indicate 
that  linear  charge  in  the  middle  part  of  the  inclusion  can 
be  successfully  calculated  by  means  of  asymptotic  equation 
(2.2.3)  for  smooth  density;  these  results  correlate  well  with 
the  classic  integral  equation  (2.2.4).  Calculations  show  that 
the  boundary  layer  makes  a  significant  contribution  near  the 
ends,  and  this  fast  change  of  density  cannot  be  governed  by 
the  asymptotic  equation  for  smooth  density  (2.2.3)  or  the 
equation  of  Pocklington.  Results  show  that  the  boundary 
layer  solution  allows  one  to  estimate  high  gradients  of  linear 
charge  near  the  ends  of  the  conductor.  Effective  zone  of  the 
boundary  layer  is  about  two  diameters  of  the  conductor. 
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Fig.  5.  Charge  density  calculated  from  equation  of  Pocklington  at  various 
aspect  ratio  (a!L-  0.05;  0.01;  0.001). 


CL 


Fig.  7.  Comparison  between  the  linear  charge  calculated  from  equation  of 
Pocklington  (dotted  line),  Eq.  (2.2.3)  (solid  line),  and  Eqs.  (2.2.3)  and  (2.2.5) 
(dashed  line)  at  a/L  =  0.05. 
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Fig.  8.  The  smooth  linear  charge  at  various  aspect  radii  (< a/L  =  0.00 1 ;  0.0001) 
for  elliptical  cross-section  b!a  =  2. 


tau 

Fig.  10.  Edge  effect  (solid  line  for  smooth  linear  charge,  dashed  line  for 
linear  with  boundary  layer);  calculations  were  conducted  for  a/L  =  0.05  and 
b!a  =  2. 


3.2.  Elliptic  cylinder 


Equations  for  the  smooth  part  of  the  linear  charge  and 
boundary  layer  were  solved  numerically,  and  dimensionless 
linear  charge  P{z/L)  (containing  smooth  part  and  boundary 
layer): 


(3.2.1) 


and  angular  distribution  lE((p)  of  charge  over  the  cross-section 
boundary  contour  were  calculated.  Fig.  8  presents  the  smooth 
linear  charge  at  various  aspect  ratios  (a/L  =  0.001;  0.0001) 
for  elliptical  cross-section  b/a  =  2.  Fig.  9  presents  the  angu¬ 
lar  distribution  ^(cp)  for  ellipses  with  ratio  b/a  =  2,  3,  and  1 
(circle).  It  should  be  noted  that  the  classical  approximation 
of  Pocklington  does  not  consider  noncircular  cross-sections. 
The  edge  effect  is  shown  in  Fig.  10  (solid  line  for  smooth 
linear  charge,  dashed  line  for  linear  with  boundary  layer);  cal¬ 
culations  were  conducted  for  a/L  =  0.05  and  b/a  =  2.  Results 
show  that  a  high  gradient  of  linear  charge  near  the  edge  can¬ 
not  be  estimated  without  calculation  of  the  boundary  layer 
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Fig.  9.  The  angular  distribution  xP((p)  for  ellipses  with  ratio  b/a  =  2,  3,  and 
1  (circle). 


term.  The  effective  zone  of  the  boundary  layer  is  about  two 
to  three  diameters  of  the  conductor. 

It  should  be  noted  that  the  2D  approach  allows  one  to 
calculate  the  angular  distribution  of  the  electric  charge  over 
the  cross-section  boundary,  but  does  not  allow  one  to  find 
the  net  charge  in  the  cross-section  z  =  const.  Finear  charge 
can  be  found  from  a  3D  formulation,  which  can  be  asymp¬ 
totically  transformed  to  an  integral  equation  over  the  axis  of 
the  conductor.  A  similar  situation  takes  place  in  the  analysis 
of  thin  films/thin- walled  inclusions.  In  some  specific  cases  a 
ID  approach  to  thin  film/inclusion  can  meet  principal  trouble 
related  to  finiteness  of  inclusion  in  the  plane  or  other  rea¬ 
sons.  This  incorrectness  of  the  boundary  value  problem  of 
the  reduced  dimension  sometimes  appears  to  be  an  artifact. 
Correct  formulation  of  the  boundary  value  problem  can  be 
obtained  on  the  basis  of  a  3D  approach.  3D  formulation  can 
take  into  account  the  finite  size  of  the  inclusion  in  the  plane, 
finite  curvature  of  the  median  surface,  non-uniform  thickness 
effects,  and  the  presence  of  additional  scale  parameter  (with 
dimension  of  length)  associated  with  the  internal  structure 
of  the  film.  All  the  above-mentioned  reasons  make  the  ID 
formulation  insufficient,  and  requires  a  more  detailed  study. 
Reduction  of  dimension  in  the  corresponding  boundary  value 
problem  is  still  possible,  but  has  to  be  implemented  more 
thoroughly.  Mathematical  methods  of  this  kind  of  reduction 
have  already  been  developed  in  solid  mechanics/acoustics. 
We  assume  that  this  advanced  asymptotic  technique  can  be 
successfully  applied  to  the  new  thin  film  effects  discovered 
in  [5]. 


4.  Conclusion 

The  electrostatic  problem  of  a  single  cylindrical  conduct¬ 
ing  inclusion  for  a  SOFC  catalyst  layers  has  been  solved  using 
the  asymptotic  decomposition  for  a  single  rod-like  conductor 
in  a  dielectric  medium.  It  was  shown  that  the  surface  distri- 
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bution  of  the  charge  can  be  described  by  the  integral  equation 
over  the  mid-line  axis  for  an  unknown  linear  charge  and  that 
this  representation  is  valid  at  the  points  far  from  the  ends  of 
the  conductor.  In  the  case  of  straight  or  a  curvilinear  inclu¬ 
sion,  the  charge  distribution  near  the  ends  is  governed  by  a 
boundary  layer  equation.  This  effective  technique  is  suffi¬ 
cient  in  describing  a  complex  system  of  various  inclusions, 
such  as  gas  diffusion  or  electrode  catalyst  layers  in  PEMFCs, 
and  a  large  number  of  non-identical  conductors  having  dif¬ 
ferent  lengths,  mid-line  contours,  and  cross-section  shapes, 
which  is  typical  for  ceramic  electrically  conductive  porous 
materials. 

As  the  next  step,  this  approach  will  be  incorporated  in 
the  complex  model  of  heterogeneous  nanostructured  mate¬ 
rial  containing  a  regular/irregular  system  of  inclusions,  which 
will  allow  the  development  of  novel  nanostructured  materials. 
It  has  been  shown,  that  derived  asymptotic  relations  are  sim¬ 
ple  enough  and  allow  extension  to  heterogeneous  medium 
associated  with  the  system  of  inclusions,  furthermore,  the 
required  level  of  accuracy  is  achieved.  It  should  be  empha¬ 
sized,  that  the  proposed  novel  approach  results  in  prediction 
of  the  properties  of  heterogeneous  material  directly,  without 
additional  simplifying  physical  hypotheses. 
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